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Abstract — This is the first of two papers describing the process of fitting ex- 
perimental data under interval uncertainty. Probably the most often encountered 
use of global optimization methods is finding the so called best fitted val- 
ues of various parameters, as well as their uncertainties, based on experimental 
data. Here I present the methodology, designed from the very beginning as an 
interval-oriented tool, meant to replace to the large extent the famous Least 
Squares (LSQ) and other slightly less popular methods. Contrary to its clas- 
sical counterparts, the presented method does not require any poorly justified 
prior assumptions, like smallness of experimental uncertainties or their normal 
(Gaussian) distribution. Using interval approach, we are able to fit rigorously and 
reliably not only the simple functional dependencies, with no extra effort when 
both variables are uncertain, but also the cases when the constitutive equation 
exists in implicit rather than explicit functional form. The magic word and 
a key to success of interval approach appears the Hausdorff distance. 

1 Introduction 

Handling experimental uncertainties lies at the heart of physics and other sciences. The 
situation is rather clear during simple measurements, for example when determining 
the mass, m, of an object under investigations. In this case the uncertainty of a single 
measurement is equal to one half of the smallest division (A) of the measuring device, 
here the balance. Consequently, we can be sure that the unknown mass certainly belongs 
to the interval [m — A/2, m A- A/2], where to is a direct readout from the instrument's 
scale. We can try to increase our confidence into the result by repeating the measurement 
several times or by using several measuring devices, maybe delivering the results with 
different precision A. By doing so, we necessarily switch to the probabilistic way of 
presenting measurement's results, namely as a pair of two numbers: the most likely value, 
TO, and its standard deviation, (t(to), expressing the uncertainty of a measurement. 

Nowadays more and more measurements belong to a class called indirect measure- 
ments, where we deduce the numerical values of various interesting parameters without 
comparing them directly to the appropriate standard. When determining the resistance 
of a sample we usually take a series of measurements, obtaining many pairs (voltage, 
current). Then, applying the Ohm's law relating those two physical quantities, we find 
the single proportionality constant, called resistance, R = voltage/current (or its inverse: 



conductance). Using many measurements, not just one, we switch naturally to the prob- 
abilistic form of result's presentation: as the most likely value and its standard deviation 
(here: R and cr{R)). 

This situation is highly unsatisfactory. First, we can never be 100% sure whether the 
partial uncertainties are indeed small, thus justifying the application of various 'error 
propagation laws'. Secondly, neither the probability density function of measurements nor 
of the result need not to be normal (Gaussian) , thus invalidating the frequent claim that 
the probability of the true value belonging to the interval [m — a, m + a] is approximately 
equal to 67%. In reality, the overwhelming majority of contemporary measuring devices 
deliver discrete, digital results. Such results obviously cannot be normally distributed. 

In this paper I present other ways of experimental data processing, using interval 
methods. Here we operate on guaranteed rather than probabilistic quantities, thus the 
results obtained on this way should also be guaranteed. In practice, the results pro- 
duced by straightforward interval implementation of the classical methods often appear 
disappointing. In this paper I point to the possible reasons and show a remedy. 

2 Bcisics of interval calculus 

From now on we will use intervals every time our numbers are uncertain. Specifically, by 
interval x we will understand a (sub) set of real numbers: 

x=\x,x] = {R^ 3 x: x^x ^x} . (1) 

The set of all intervals is usually written as IR. The familiar real numbers can be identified 
with intervals of type x = [x,x], i.e. having their lower and upper bound equal to each 
other, and called degenerate intervals, thin intervals, or singletons. This way the set of 
intervals may be regarded to be an extension of real numbers: C IK. 

It is possible to define the arithmetic operations on intervals: 

addition: x + y = [x + y, x + y\; 
subtraction: x — y = [x — y, x — yj; 

multiplication: x • y = [minZ, maxZ], where Z = \^xy, xy, xy, xy}; 
division: x/y = x • i, where i = [l/y, l/y\ , ^ y. 

One can understand the right hand sides of the above definitions as intervals certainly 
containing the results of the corresponding operations when the first operand is arbitrarily 
chosen from interval x and the second from interval y. One can also easily check that 
we recover precisely the ordinary results whenever x = x and y = y.lt has to be stressed, 
however, that in computer calculations it is necessary to properly round every result. The 
so called outward rounding is appropriate here, that is the lower bound (x) should be 
always rounded towards — cxd and the upper bound (x) should be always rounded towards 
+00. Only then we can be sure that our results are certain (guaranteed), i.e. containing 
the true result with probability 1. A warning is in place here: never treat the interval as 
a random variable uniformly distributed in its interior. If this was so, then the sum of 
two intervals would have to represent a triangularly distributed random variable, what 
is - obviously - not the case. 



Since intervals are sets, then it is natural to operate on them as on sets, not only as on 
numbers. Consequently we can compute the intersection of two intervals (it may happen 
to be an empty set, conveniently written as [+00, —00] - an illegal interval), or their set- 
theoretical sum (union). The union of two non-empty intervals is never empty but not 
necessarily is an interval again. This may be sometimes troublesome, so it is desirable to 
introduce the notion of convex hull, or shortly hull of two intervals as a smallest interval 
containing them both: 



It is remarkable that the finite number of arithmetic operations (performed on inter- 
vals' endpoints) is all we need to process uncountably infinite subsets of real numbers. 
The nice feature of interval computations is its straightforward extensibility to operate 
on vectors, matrices, or real-valued functions as well. 

The essence of interval calculations is to deliver guaranteed results, but sometimes we 
obtain intervals which are wider^ than necessary. The overestimation is usually hard to 
asses (and may be quite large!) but will happen almost surely when in a given expression 
any variable appears more than once; for example: [0, a] — [0, a] = [—a, a] and not = [0,0] 
— as one might expect. This phenomenon is called a dependency problem. In particular, 
we have: 



So, the order of calculations matters. Fortunately, the overestimation usually vanishes 
when the widths of interval operands go to zero. Exceptions do exist but it is not easy 
to present a practical example. 

The vectors with interval components are usually called boxes and it is easy to see why. 
For more information on interval computations the reader is referred to other sources [1] . 

2.1 Common features of interval algorithms 

Interval-oriented algorithms usually operate on lists of boxes and treat them accord- 
ingly to their current status. Good boxes are retained for further reference, while bad 
boxes are discarded as soon as possible. The third category, called pending, unresolved 
or undetermined, is the most interesting and as such is the main object of processing. 
The general idea is to start from a single big box, suspected to contain the solution (or 
solutions) and therefore initially labelled as pending. Pending boxes, one by one, are 
divided into smaller parts (usually two) and disappear from the list. Offspring boxes are 
subjected to one or more tests aiming at determining their current status. Some of them 
are recycled back onto pending list, while bad are immediately discarded, and good are 
collected separately as a part of the final result. The procedure terminates when either 
the list of pending boxes is empty or contains only 'small' boxes. On exit, the result is 
a union of good boxes, possibly appended with pending ones. Empty output is a proof 
that the solution(s) of our problem, if any, are located outside the initial box. 

Of course, the exact meaning of bad and good depends on context. 



The width of interval is a positive real number: w(x) = x — x. 



hull(x,y) = D(x,y) = [mm{x,y), max(a;,?;)] 



(2) 



X (y + z) C xy -|- xz 



(3) 



3 The problem of experimental data fitting and its solution 



Suppose we have a theory T, characterized by k unknown parameters pi, . . . , pfc, and A'' 
{N > k) results of measurements mi , . . . , rriN , each taken at different values of some 
well controlled variables, for example the varying temperature or magnetic field. Let Xj 
denotes the value(s) of controlled variable(s) (environment) during the measurement nij. 
In what follows we will use the simplified notation: 

r(pi, . . . ,pfc;xj) = tj, j = l,...,N 

where the upper line describes the theoretical outcomes of the experiment and the 
lower one - actual experimental results. Usually each theoretical outcome tj is crisp^ 
when all the arguments of T (including x^) are crisp. Contrary, the measurements rrij are 
always uncertain and will be considered from now on to be intervals (or more generally: 
interval vectors). 

The task for an experimentalist is to adjust the values of all unknown parameters 
Pi, . . . ,Pfc, in such a way that every measurement nij differs as little as possible from 
the corresponding theoretical prediction tj. In fact, we would be happy to finish the 
procedure observing 

t,nm,^0, Vj e {l,2,...,iV} (5) 

(in simple words: the theoretical curve passes through all the experimental points) and 
with unknown parameters, p's, as narrow as possible. 

Note that the requirement tj C nij for every j, although tempting, is going too far: 
our results (if they exist) would be severely underestimated. Yes, we want to be accurate, 
but we can't afford to tolerate underestimates. 

It is out of scope of this article to elaborate the details of many existing particular 
procedures [2]-|14j aimed to maximally contract the initial box of unknown parameters 
Pi X • • ■ X Pfc in such a way that all the relations (O are satisfied. Some of those methods 
(see [9]) deliver a single box being an interval hull of solutions, while others, more time- 
consuming, output more boxes, covering with certainty the solution set. Technically: we 
call a box bad when at least one of inequalities ([5]) is violated at all its internal points 
and pending otherwise. In rare cases it happens to encounter tj C rrij for all j's - such 
box is good and requires no further processing. 

3.1 Advantages and disadvantages of the rigorous solution 

The practical implementations of the above described approach are still rare [TSl [H] 
and one may wonder why. The unquestionable mathematical rigor of a procedure is 
certainly among its major advantages. Unfortunately, for very practical reasons, this is 
also its weak point. And here is why: 

• the routine requires guaranteed intervals as measurements, i.e. the ones containing 
true outcomes of an experiment with probability equal exactly to 1. It is impossible 
to satisfy this requirement in laboratory practice since measurements usually have 
the form m = nio ± cr(m), where toq is a mean value and a{m) - its standard 



^We call the object crisp (point, point-wise, point-like) in contrast to the one having interval char- 
acter, be it a real number, vector or matrix. 



deviation. Nobody knows (or cares) what is the distribution of m, even its support is 
usually unknown (well, some physical quantities, like mass or absolute temperature, 
have to be non- negative). Consequently, taking intervals [toq — a, mg + a] as input 
data, we will most likely finish with empty set of results. This is because at least 
one-third of N relations ^ has to fail, for arbitrary (even crisp!) set of unknown 
parameters {pi,. ..,pk}; 

• the quick and dirty fix, coming to the mind in the above situation, is to use wider 
intervals, like [mo — Su, uiq + 3a], as 'almost guaranteed' input data. Depending 
on the individual luck, this trick may or may not work. If it doesn't then we are 
left without even a foggy idea what are the values of our unknown parameters. 
If it works then uncertainties of so obtained parameters often seem exceedingly 
large, very pessimistic, at least when compared with those obtained on other ways. 
Besides: are we entitled to 'scale back' the obtained uncertainties dividing them 
all by 3? The only honest answer is no; 

• sometimes, when the fix works, the uncertainties of unknown parameters appear 
unrealistically, not to say suspiciously, small. This will surely happen when among 
our data there is at least one element which should be labeled as a 'near outlier'. 
It may be due to the undetected data transmission/recording error, power line 
fluctuations, or whatever. 

• another approach is to require fitted curve to pass through at least Ni < N exper- 
imental points, without any prior indication which points are to be preferred. The 
sensible choice for A^i is iVi « 2N/3. We will not discuss further this idea. 

On the other hand, when all our data are credible but the solution set is empty - our 
theory T must be flawed. This could be a very strong statement, but is not. There is 
one more possibility: the experimental uncertainties are seriously underestimated, delib- 
erately or otherwise. So, experimenters, be warned: cheating will be severely punished 
by interval methods and will bring you nowhere! 

Concluding: while the rigorous approach has many advantages, we definitely need 
something else. 

3.2 Back to the basics 

The situation is especially upsetting when experimental data look fine, and we are 
sure the theory is correct, yet the interval routine returns no answer at all. What can we 



In classical data analysis, the LSQ method (Least Squares) is most commonly used 
to find unknown parameters. In short, its essence is to minimize the quantity called chi 
squared: 



where p is a (crisp) vector of k unknown parameters - arguments of the function x^- We 
are looking for such a vector p* as to have x^iP*) = "lin- It is well known that LSQ 
method never fails and always produces some results, even for completely wrong theory. 

The LSQ method is due to Carl Friedrich Gaufi (1777-1855) and was originally in- 
vented by him around 1794. Later on, in 1809, the same author gave it solid statistical 



do? 




(6) 



interpretation. We will not proceed with statistical properties of and LSQ method in 
general, instead we will rather concentrate on the original idea. And it was like that: 

Due to unavoidable experimental uncertainties (then called errors), it is unlikely to 
draw a line representing theoretical outcomes as a function of environment (x) and have 
the experimental points mj (x) lying precisely on this curve at the same time. Note that 
the notion of standard deviation of the measurements was practically unused those days, 
so the measurements were just numbers. This means there were no denominators in 
formula Varying the unknown parameters pj's, we deform the theoretical curve so 
as to have it running as closely to the experimental points as possible, thus solving the 
problem. 

Remark: 

Gaufi could not speak about the distance (between a theoretical curve and 
experimental points), at least not in a strict mathematical sense, since the 
notion of distance was introduced into mathematics many years later, by 
another German mathematician, Felix Hausdorff (1868-1942), and, indepen- 
dently, in 1905, by Romanian mathematician Dimitrie Pompeiu (1873-1954) 

m- 

Having known the notion of distance, Gaufi would almost certainly use it in his early 
idea of LSQ method. Indeed, the formula ([6|) is nothing else as squared Euclidean distance 
in A'^-dimensional space, where aj serves as a unit length in direction j. 

4 Interval version of LSQ method 

The problem of global optimization is well studied. Interval researchers have also con- 
tributed their share to this field. They have already devised many interval-thinking- 
inspired, rigorous procedures aiming to solve such tasks. Minimizing is obviously one 
of them. Simple optimization cases, in LSQ sense, were investigated since at least 1990 
|18j . Several closely related problems have been successfully solved using rigorous interval 
algorithms, yet the specific difficulties, recalled in Sec. 13.11 still have not been addressed 
satisfactorily. 

The enormous potential of interval methods for uncertain data processing has been 
recognized some time ago (1993) by Walster and Kreinovich [T^. Kosheleva and Kreinovich 
(1999) pointed that the cost of interval approach only slightly exceeds the traditional 
(probabilistic) one [20J. Muhoz and Kearfott (2001) have shown that non-smooth cases 
should be essentially no more difficult than regular ones [TJ. Only one year later (2002), 
Yang and Kearfott heralded 'new paradigm' and 'new way of thinking about data fitting' 
|21j . Unfortunately, they did not show how to practically implement their ideas other 
than for linear equation set. Finally, (2005) Zhilin [T3j found the solution for the case 
when measurements are (multi)linearly related to the environment. 

4.1 Correct interval form of 

How to 'translate' the fundamental formula ® into its interval counterpart? The first 
idea is to simply replace real-valued measurements mj by their interval representations 
m^ , apply the same trick to theoretical predictions tj tj , and retain real- valued aj 's. 



Consequently, the functional becomes interval-valued as well. So far, so good, but 
wait. 

In the light of what was said before, we have to rewrite in terms of distance 
between experimental and theoretical (predicted, simulated) points. Yet the expression 
nij — tj, nor even \vi\j — tj|, is not the correct mathematical distance! There are two 
reasons for that: 

• it is not a real number, and 

• the implication (m = t) (m — t = ) is false. 

This naturally raises the question: what is the distance between two intervals, here 
m and t? Following Hausdorff, we will use his metrics, adapted for IK space by Moore 
as: 

d (a, b) = max (| a — 6 | , | a — 6 | ) . (7) 

Before continuing, let us only note that another function d! (a, b) = C • d (a, b) , where 
C > is a fixed number, is a correct distance, too. Additionally, the mathematical 
distance is not expressed neither in miles nor millimeters, or in any other units - it 
is simply dimensionless. This is important, since using the strict mathematical distance 
between intervals we will be able to fit our experimental data obtained from two (or more) 
completely different experiments, provided the sets of unknown parameters, relevant to 
each experiment separately, have non-empty intersection. Simply speaking - it is possible 
to fit several curves simultaneously, in a single run. 

Remarks: 

The Euclidean distance between theoretical and experimental points is not 
the only thinkable distance. We will discuss other iV-dimensional metrics in 
the other article. We will stick, however, to the Moore-Hausdorff distance as 
one-dimensional metrics in IR, no matter that there are other choices, see for 
example [13] and [H] . Note also that Moore-Hausdorff distance specialized to 
the real line simplifies to the familiar form: d(a, b) ~ | a — & |, as one might 
expect. 

4.2 Analyzing term by term 

Consider now again a single term in What is the distance between the predicted 
value t and the true value m* {m* is a real number, not interval)? All we know about m* 
is that it satisfies the double inequality: m ^ m* ^ m, with exact value of m* remaining 
unknown. According to the definition ^ we have {m* — [to*, to*] — m*): 

d(t, TO*) =max(|i-TO*|, |Z-TO*|) e R (8) 

We don't know which internal point of m is equal to to*. Suppose for a moment that 
to' coincides'^ with one of the endpoints of m, say to' = to. Call the current distance 
d(t, to' = to) = ^. What happens when to' gradually moves to to - the other endpoint 
of m? 

If t and m are disjoint then ^ will linearly increase or decrease, depending on the 
relative position of t and m on a real axis. We will finally get either d (t, to) = -I- w(m) 



•^For a moment we will be dealing with some m', not even necessarily contained in [m, m], rather 
than with m*. This is just for purity: m* is a fixed number and we want m' to be variable. 



or d (t, m) = w(m) (this must be a positive number as t n m = 0). Consequently we 
have bounded the true distance d{t,m*) with uncertainty |d(t,TO) — d(t,rn)\ = w(m). 
This is indeed very remarkable result, especially when compared with the so called natural 
interval extension of a basic building block of namely the expression 't — m'. There 
we always have w(t — m) = w(t) + w(m) > w(m). In our approach the uncertainty of 
a distance in question never exceeds the uncertainty of an individual measurement and 
- surprisingly - it does not depend on current accuracy of unknown parameters. 

It remains to show how this result changes when t and m overlap. Previously, thanks 
to the disjointness of t and m, only one argument of ^ was 'active' at all times, meaning 
that its value defined the distance. Now, at some point m' (not necessarily m*), the two 
arguments can exchange their roles and the other one may become 'active'. If this happens 
(it may not) then the direction of change of ^ will change too, thus 'un-doing' the already 
acquired change. Therefore the final change of ^ cannot even reach w(m). 



5 Main result 



The analysis just presented is intuitive and easy to follow, but still crude. It says nothing 
about the values of bounds, being limited only to their separation. This is not enough 
to be useful in practice. In particular, it is easy to show that m' at which the two 
arguments of ^ are equal to each other is m' = (t + t) /2 - the center of interval t (but 
only when m' also belongs to m, otherwise there is no such point and switching of roles 
does not occur). For such m', the distance d (t, m') — (t — i)/2. With this fact in mind, 
we are able to derive the exact bounds for the distance p, between a predicted value t of 
an experimental outcome, under known environment x, and the true value m* G m, in 
the same circumstances x. They are following: 



when c(t) G m; 



when c(t) ^ m: 



lower bound: p — ^w(t) 

upper bound: p = max [d{t,rn), d(t,m)] 

lower bound: p = min [d (t, m), d (t, m) ] 

upper bound: p — max [c? (t, tti), c? (t, to) ], 



(9) 



(10) 



where c(-) stands for the center of its interval argument, here c(t) — ^ {t + tj, and d (•, •) 
is a Moore-Hausdorff distance between intervals. Please note that generally 



d (t, m) 7^ max [d (t, to), d (t, to) ] , 



(11) 



so the above bounds cannot be written in a more compact form. 

In addition, we have introduced a new symbol, p, for the mathematically correct 
distance between the true measured quantity and its interval theoretical estimate, p 
stands for its interval enclosure. The intuitively shown relation w(p) ^ w(m) remains 
true. 



6 Discussion 



In conclusion, the recommended form of interval version of ^ functional is: 

p(tj-(pi,...,pfc), mj) j^ 

where p = is given by jl]) and (fTO|) . respectively. To recover the unknown pa- 

rameters pi, . . . ,pfc, one has to find a global minimum of (|12p with respect to those 
parameters. This task may be accomplished using procedures developed mostly by Luc 
Jaulin and Eric Walter (notably SIVIA algorithm, set-inversion methods) as well as ideas 
first put forward by Shary PS]. No matter which approach will be adopted, it is clear 
that no crisp values can be ever obtained, just because all m's are intervals. On the other 
hand, this very feature is highly desirable, since the uncertainties of searched parameters 
are evaluated very credibly and as nearly a side effect, with no extra effort. This obser- 
vation sheds new light not only on the problem of experimental data fitting in general, 
but also substantially changes our perspective on reliable estimates of uncertainties in 
indirect measurements. 

Once the measurements are completed, the values of all m's, as well as their widths, 
are fixed. This makes possible to use widths of measurements as the natural unit lengths 
in every direction of N possible. That is why m's are present in denominators of all 
components of the sum (fT^ . When we speak of distances this is the only choice of appro- 
priating individual weights to all measurements. There is no space left for arbitrariness, 
as it sometimes happens in other versions of the so called weighted LSQ regression. 

As the measurements are fixed during calculations, the interval enclosure p of the 
distances between the true unknown values m* and their theoretical predictions tj, is 
perfect and cannot be further improved. Moreover, those widths are uniquely determined 
by the accuracy of the real measurements, not guessed or necessarily subjectively evalu- 
ated by a human expert. The optimal widths are very fortunate, since the lower is the 
width of the interval extension of a function being minimized, the more precisely the 
global minimum may be located. At least at this respect our procedure definitely beats 
the natural approach. 

In the second part of this work a connection between interval hulls of solutions and 
statistical description of their uncertainties will be demonstrated and discussed. 

Acknowledgments 

This work was done as a part of author's statutory activities at the Institute of Physics, 
Polish Academy of Sciences. 

Bibliography 

[1] The website http: / / cs.utep.edu/interval-comp / main. html' is an excellent and up-to- 
date entry point to the wonderful world of interval computations. 



[2] Luc Jaulin and Eric Walter, Guaranteed Nonlinear Parameter Estimation via Inter- 
val Computations, Conf. on Numerical Analysis with Automatic Result Verification, 
Lafayette, Feb. 25th-Marcli 3rd, 1993, pp. 61-75 

[3] Luc Jaulin and Eric Walter, Guaranteed nonlinear parameter estimation from 
bounded-error data via interval analysis, Math, and Comput. in Simulation 35, 
123-137, 1993 

[4] Luc Jaulin and Eric Walter, Set Inversion via Interval Analysis for Nonlinear 
Bounded-error Estimation, Automatica 29, 1053, 1993 

[5] L. Jaulin and E. Walter, Guaranteed Parameter Bounding for Nonlinear Models with 
Uncertain Experimental Factors, Automatica 35, 849-856, 1999 

[6] Luc Jaulin, Interval constraint propagation with application to bounded-error esti- 
mation, Automatica 36, 1547, 2000 

[7] Humberto Munoz and R.B. Kearfott, Interval Robustness in Non- 
smooth Nonlinear Parameter Estimation, unpublished preprint, 
http://interval.louisiana.edu/preprints/2001_robustness.pdf 

[8] L. Jaulin and E. Walter, Nonlinear Bounded-Error Parameter Estimation Using 
Interval Computation, in: Granular computing: an emerging paradigm, Physica- 
Verlag GmbH Heidelberg, pp. 58-71, 2001 

[9] Marek W. Gutowski, Prosta dostatecznie gruba, Post§py Fizyki, 53(4), 181 — 192, 
2002 (in Polish) {Fat enough straight line. Advances in Physics, bimonthly of Polish 
Physical Society) 

[10] Voschinin Alexander, Tyurin Alexander, Interval identification of time series pa- 
rameters using readings with bounded errors, 5th International Scientific- Technical 
Conf. ProcessControl, Pordubice 2002, paper R-210, 7 pages 

[11] M.H. van Emden, Using the duality principle to improve lower bounds for the global 
minimum in nonconvex optimization. Second COCOS workshop on intervals and op- 
timization, 2003 (published version: Using Propagation for Solving Complex Arith- 
metic Constraints, http://arxiv.org/abs/cs/0309018) 

[12] I. Braems, N. Ramdani, A. Boudenne, L. Jaulin, L. Ibos, E. Walter, and Y. Can- 
dau. New set- member ship techniques for parameter estimation in presence of model 
uncertainty, Proc. of the 5th Int. Conf. on Inverse Problems in Engineering: Theory 
and Practice, Cambridge, UK, 11-15 July 2005 

[13] Sergei I. Zhilin, On Fitting Empirical Data under Interval Error, Reliable Comput- 
ing 11, 433-442, 2005 

[14] Maarten van Emden, Constraint-Driven Global Optimization, 13th International 
Symposium on Scientific Computing Computer Arithmetic and Verified Numerical 
Computations SCAN'2008, El Paso, Texas, September 29 - October 3, pp. 144-145, 
2008 

[15] L. Jaulin, J-L. Godet, E. Walter, A. EUiasmine, and Y. Le Duff, Light scattering 
data analysis via set inversion, J. Phys. A: Math. Gen. 30, 7733-7738, 1993 

[16] Maelenn Aufray, Adrien Brochier, and Wulff Possart, Set Inversion via Interval 
Analysis applied to dielectric spectroscopy , the talk given at SWIM 2008, June 19- 
20th, Montpelier, France 



[17] T. Barsan, D. Tiba, One hundred years since the introduction of the set distance 
by Dimitrie Pompeiu, in: System modeling and optimization, Springer, New York, 
pp. 35-39, 2006 

[18] Svetoslav M. Markov, Least-square approximations under interval input data, Con- 
tributions to Computer Arithmetic and Self- Validating Numerical Methods, C. Ul- 
rich (editor), J.C. Balzer AG, Scientific Publishing Co. (C) IMACS 1990, pp. 133-147 

[19] G. William Walster and Vladik Kreinovich, For Uknown-but- Bounded Errors, In- 
terval Estimates are Often Better Than Averaging, ACM SIGNUM Newsletter 31, 
6-19, 1996, http://www.cs.utep.edu/vladik/1993/tr93-31b.ps.gz 

[20] Olga Kosheleva and Vladik Kreinovich, Error Estimation for Indirect Measurements: 
Interval Computation Problem Is (Slightly) Harder Than a Similar Probabilistic 
Computational Problem, Reliable Computing 5, 81-95, 1999 

[21] Jie Yang and R. Baker Kearfott Interval Linear and Nonlinear Re- 
gression - New Paradigms, Implementations, and Experiments or New 
Ways of Thinking About Data Fitting, talk given at the Seventh 
SIAM Conference on Optimization, May 20-22, 2002, Toronto, Canada, 
http://interval.louisiana.edu/preprints/2002_SIAM_minisymposium.ps 

[22] R.E. Moore Interval Analysis Prentice Hall, Englewood Cliffs, NJ, 1966 

[23] Francesco Palumbo and Antonio Irpino, Multidimensional Interval- 
Data: Metrics and Factorial Analysis, Applied Stochastic Models 
and Data Analysis, Brest, France May 17- 20, 2005, pp. 689-698, 
http://webhouse.unimc.it/economia/repo/39/689.pdf 

[24] Antonio Irpino and Rosanna Verde, Dynamic clustering of interval data using a 
Wasserstein-based distance. Pattern Recognition Letters 29, 1648-1658, 2008 

[25] Sergei P. Shary, A Surprising Approach in Interval Global Optimization, Reliable 
Computing 7, 497-505, 2001 



